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EMPIRICAL MODEL OF GLOBAL THERMOSPHERIC TEMPERATURE 
AND COMPOSITION BASED ON DATA FROM 
THE OGO-6 QUADRUPOLE MASS SPECTROMETER 

ABSTRACT 

An empirical global model for magnetically quiet conditions has been 
derived from longitudinally averaged N , O, and He densities by means of 
an expansion in spherical harmonics. The data were obtained by the 
OGO-6 neutral mass spectrometer and cover 'the altitude range 400 to 

*4 

600 km for the period 27 June 1969 to 13 May 1971. The accuracy of the 
analytical description is of the order of the experimental error for He and 
O and about three times experimental error for N 2 , thus providing a reason- 
able overall representation of the satellite observations. Two model 
schemes are used: one representing densities extrapolated to 450 km and 
one representing densities extrapolated to 120 km with exospheric temperatures 
inferred from N 2 densities. Using the best fit model parameters the global 
thermospheric structure is presented in the form of a number of contour plots. 
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EMPIRICAL MODEL OF GLOBAL THERMOSPHERIC TEMPERATURE 


AND COMPOSITION BASED ON DATA FROM 
THE OGO-6 QUADRUPOLE MASS SPECTROMETER 

I. INTRODUCTION 

This report presents a global thermosphere model for quiet magnetic con- 
ditions based on in situ measurements of N 2 , O, and He obtained with the neutral 
particle quadrupole mass spectrometer [Carignan and Pinkus, 1968] carried 
aboard the OGO-6 satellite. The satellite was launched June 5, 1969 into an 82° 
inclination orbit with a 398 km perigee and 1100 km apogee and was operated until 
July, 1971. The reduction of the raw data to ambient densities makes use of the usual 
thermal transpiration equation in a moving coordinate system [Schultz et al., 1948; 
Horowitz and LaGow, 1957] and further includes corrections to the N 2 densities 
as a result of CO contributions to the mass 28 peak, and corrections to atomic 
oxygen densities as a result of oxygen surface adsorption, recombination, and 
desorption [Hedin et al., 1972a] . For each gas species a density value is 
determined every 9 seconds. 

The uncertainty in N 2 and O densities as a result of noise in the data is 
always about 2% to 4%, but for He the uncertainty varies greatly (depending on 
the exact conditions) because of the generally low signal to noise ratio. Near 
perigee the uncertainty in density due to possible systematic errors in back- 
ground subtraction and gas-surface interaction corrections is generally about 
2% for N 2 and 6% for O (although it can be several times larger under very low 
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density conditions). These systematic errors increase with altitude, particularly 
above 500 km. No densities are used for which the total uncertainty (random 
and systematic) is estimated to exceed 25% for N 2 and O and 50% for He. In 
addition there is a laboratory calibration uncertainty of 10% to 15%. 

The model described here is a more recent version of a model presented 
byHedinetal. [1972b] at COSPAR. The principal changes are inclusion of: 8 
hour local time components in the model formula; data from days with slightly 
higher magnetic activity; and data from the second year of satellite operation. 
The model is intended to provide in concise form as accurate a representation of 
the thousands of measured ambient densities as possible. At the same time, the 
model reveals many features, both familiar and novel, in the global distribution 
of the various gases which are inherent in the data but would not otherwise be 
easily seen. 
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II. DATA SELECTION AND COVERAGE 

To eliminate strong magnetic activity effects data are selected 
from days which have a daily magnetic index, A p , of 7 or less and a 3 hour 
magnetic index, a p , less than 12 for that day and 6 hours earlier. The quiet days 
which have at least 8 orbits of good data reasonably distributed throughout the 
day are listed in Table 1. All the data obtained on a given day are grouped into 
5° latitude bands (distinguishing also between data taken on north or south bound 
passes since these represent drastically different local times) and averaged over 
a 24 hour (universal time) period. All averages used in the model are based 
upon at least 20 density points. The averaging procedure eliminates 
longitudinal/UT effects [Hedin and Reber, 1972] in the data and reduces the 
amount of data handled, but does not eliminate local time effects since the orbit 
plane changes quite slowly with respect to the sun (2°/day). 

For He there is a special data selection problem because it is observed 
that data from the second year of operation frequently differ by large factors 
from a model using first year data only. Raw data obtained during these periods 
are observed to fluctuate erratically in magnitude and the peak shape is not 
normal. Thus, an iterative procedure is used in fitting the data in which density 
points which differ from the model by more than a factor of two are eliminated 
and the remaining points refit. 

The latitude and local time coverage of the data can be seen in Figures 1 
and 2. In one year the latitude of perigee made three cycles around the globe 
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and (effectively) five diurnal cycles. Values of the mean (taken over three solar 
rotations) FlO.7 flux vary from 108 to 168 with an overall average near 150. 
Departures of the daily F 10 7 flux from the mean vary from -40 to +60. The 
overall average Ap value is 4. The number of days with good data available in the 
second year of operation is much less than in the first year because of telemetry 
and spacecraft problems. 
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HI. MODEL FORMULA AND DATA FITTING 


A. 450 km Densities 

The composition is described in terms of density variations at 450 km, with 
an altitude dependence which assumes isothermal conditions (above 400 km) at 
any given time and location. Using a least squares technique, the average density 
data are fit with the following formula based on a spherical harmonic expansion 
(in geographic latitude -local time coordinates): 


n (z, L) = n (450) exp [G (L) - 1] exp 


Mg (z - 450) (R e + 450) 
R T a) (L) (R e + z) 


where 


( 1 ) 


n = number density of a particular gas species [cm - 3] 
z = altitude (km) 

L = variable representing the geographical and geophysical parameters 
on which it is assumed density depends 
M = molecular weight (gm/mole) 

R e = radius of the earth (6356.77 km) 

g = acceleration of gravity (8.5529 x 10 ” 3 km/sec 2 at 450 km) 

R = gas constant (8.314 x 10 -3 gm-km 2 /mole-sec 2 -deg) 

Tod = exospheric temperature (°K) 

G = spherical harmonic expansion as described in text 
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The function G (L) contains the basic information on the variation of density with 

selected geographic and solar activity parameters and is given in general as: 

00 

G(L) = 1 +2] a n (td- F 10.7' V P n ^ 

n= 1 

co n 

+ ^ ^ td ’ Fl0 -7’ V P n ^ C ° S " T ) 

n= 1 m= 1 

oo n 

+ ^~! 22 b n ^d- P l 0 , 7 ’ A p) P n ( 0S > sin (m O) T) 
n= 1 m s 1 

where 

a™ = parameters of the spherical harmonic expansion 
which are assumed to vary with day of the year 
(t d ), solar flux (F 10 7 ) and magnetic activity (A p ) 

P™ = Legendre polynomials (not normalized) 

6 = colatitude 


co = angular rate of earth rotation 
T = local solar time 

The merit of an empirical model based on a spherical harmonic analysis 
is that these functions form a complete set which can, in principle, represent 
any degree of complexity in the data by systematically increasing the number of 
coefficients used. Furthermore, the spherical harmonics are approximate 
eigenfunctions in the thermosphere so that relatively few terms should be re- 
quired, and analysis in terms of these components may aid in understanding the 
origin of various thermospheric variations [ Reber, 1971; Mayr and Volland, 1971; 
Volland and Mayr, 1972; Reber and Hays, 1972] . 
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The specific expansion of Equation 2 used to fit the data for the current model 


is as follows: 


G = 1 


(time independent) + a 2 o P 20 + a 40 P 40 

(primary F 10 ? ef feet) + F Q 

(magnetic activity) + (k® Q + k * 0 P 2Q ) AA 

(symmetrical annual) + (c * 0 + c* Q P 2Q ) cos f) (t rf - tgj) 

(symmetrical semiannual) + (c ^ 0 + e| Q P 2Q ) cos 2 fi (t d - t“) 

(asymmetrical annual) + (cj Q P 1Q + c* Q P 30 + cj fl P s0 )(l + Fj) cos G(t d - t^) 

(asymmet rical semiannual)+Cj Q P 1Q cos 2 £1 (t rf - t^j*) 

(diurnal) +fa P +a P +a P 

v ' r 1 11 r ll + 31 31 + 51 r 51 

+ (cjj P n + P 21 ) cos n (t d - t^ 1 )] (1 + Fj) cos cot 

+ ^11 P 11 + ^31 P 3 1 + ^51 P 51 

+ ( d li P 1 1 + d 2i p 2i) cos n (td - t? 5 )] (! + Fj) sinwr 


(semidiurnal) + [a 22 P 22 + c‘ 2 P 32 cos Q (t d - t«)] (1 + F,) cos 2a ,r 

+ [b 22 P 22 + d 32 P 32 cos fi (t d - t^)] (1 + Fj) sin 2 a it 
( terdiurnal) + a 33 ? 33 (1 + Fj) cos 3 co r 

+ b 33 P 3.3 ( 1 + F i) sin 3 WT ( 3 ) 
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where 


P = P m (cos -0 ) 

nm n v ' 

F o = f oo AF+f ooCAF ) 2 + f^zrF 
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= F 

- F 

10.7 

1 10. 7 

“ p 
r 10 . 7 

- 150. 


A A = A p - 4. 

F 10 ? = 10.7 cm flux one day earlier (10~ 22 watt - m' 2 - cps’ 1 ) 

F 10 ? = 3 solar rotation average of F ]Q 

Ap = 24 hour average of magnetic index, a p , with a six hou r 1 ag 
0 = 2 tt/ 365 (days' 1 ) 
co - 2 n/24 (hr s" 1 ) 

t d = daycount (days) 
r - local solar time (hrs) 

The significance of the various terms in Equation 3 is indicated briefly in paren- 
thesis at the left. The expansion includes local time independent, diurnal, semi- 
diurnal, and terdiurnal terms plus terms representing the annual, semiannual, 
solar flux, and magnetic activity variations. For convenience, terms higher than the 
eight hour local time harmonic are not included. However, it has been shown that the 
first three harmonics are sufficient to fit incoherent scatter data reasonably well 
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[Salah and Evans, 1972] . The form of the F 10 7 correlation used is similar to 
that of Jacchia [1971] except that allowance is made for a possible non- 
linearity in the correlation with the 27 day F 10 7 variation as suggested by in- 
coherent scatter data [Waldteufel and Cogger, 1971] . Following Jacchia it is 
also assumed that changes in the diurnal and annual temperature amplitudes 
resulting from F variations are strictly proportional to the changes in mean 
temperature, but the same assumption is not necessarily true for thermospheric 
densities which depend strongly and nonlinear ly on temperature through the 
hydrostatic equation. Thus it is necessary to introduce the parameter f when 
determining the 450 km densities for N 2 and O (fj did not differ significantly 
from unity for He). When determining exospheric temperatures and 120 km den- 
sities, as described in the next section, f l is assumed to equal one. The cor- 
relation with A p is assumed to be linear, as the range of A p included in this 
quiet time model is quite small. While the A p used in the fitting is a daily 
average based on a six hour lag as adopted in the Jacchia models, 
the difference in calculating densities using the normal daily A p is insignificant. 
The possibility of a latitude variation of the magnetic activity effect was included. 
Combining Equations 1 and 3 there are 37 unknown parameters to be determined 
by fitting the data. 

The temperature (T^,) needed in Equation 1 is determined from the N 2 data 
in a manner to be described below, but is not critical to determining 450 km 
densities since this altitude is near the average altitude of the measurements. 
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Tests were made using both temperatures from Jacchia [1965] and a constant 
temperature to verify that the 450 km density determination is relatively in- 
sensitive to the exospheric temperature. 

B. Inferred Exospheric Temperatures and 120 km Densities 
Adopting a method analogous to that employed by Jacchia [1965] , the tem- 
perature and N 2 density are fixed at 120 km and the exospheric temperature 
necessary to calculate the measured densities is then determined using analytic 
density profiles suggested by Walker [1965] based on a temperature profile 
suggested by Bates [1959] : 

T (z, L) = T* (L) - (T* (L) - T 120 ) exp [- s (z - 120)] (4) 

n (z, L) = n (120) D (z, (L), s) (5) 

/ , . \i+o + r 

D (z, T ot (L), s) = ( ) exp [- o' y £] • (6) 

\1 - a e-°V 

where 

cr = s + .00015 

a = 1 - T /T (L) 

i = (z - 120) (R + 120)/ (R + z) 

7 = M g /(or RT a (L)) 

T ro (L) = T, G (L) 

T 120 = temperature at 120 km (355°K) 
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s = temperature gradient parameter 

g e = acceleration of gravity at 120 km (9.44663 x 10“ 3 km/sec 2 ) 

= average exospheric temperature 

The temperature calculation is based on the N 2 density, rather than O or 
He, because this species is probably most nearly in diffusive equilibrium [Reber, 1971; 
Mayr and Volland, 1971; Mayr and Volland, 1972a] when dynamical processes 
are considered. In addition, variations in the eddy diffusion will have the least 
effect on the nitrogen concentration because its molecular weight is approximately 
the mean molecular weight at the turbopause. However, the average temperature 
(T ot ) needed to fit the measurements is not arbitrary, because the measurements 
are made at different altitudes, and there is an average temperature which will 
best fit the density data. The capability of determining an average temperature 
independent of the 120 km boundary conditions was maintained by allowing the 's' 
parameter to be determined by the least squares fit. Combining Equations 3, 5, 
and 6 there are 37 unknown parameters to be determined by fitting the nitrogen 
data, including T ro , s, and the parameters in function G. 

It would be possible in principle to determine independently the global 
temperature distribution from scale height information alone and a preliminary 
attempt indicated that results would be in rough agreement with the temperature 
distribution obtained using the diffusive equilibrium model of Equation 5. How- 
ever, the uncertainty limits were so large that it was doubtful that useful addi- 
tional temperature variation information could be obtained from the scale height. 



Using the temperatures inferred from the nitrogen data, the average den- 
sities and variations at 120 km needed to reproduce the measurements of O and 
He are determined by using a slight modification of Equation 5: 

n (z, L) = n (120) D (z, (L), s) exp [G (L) - 1] (7) 

where the density variations are again introduced through the G function of 
Equation 3. There are 36 unknown parameters to be determined including 
n (120). 

It should be noted that the 450 km density model and the diffusive 
equilibrium model starting at 120 km can (in principle) represent the measured 
densities equally well, as the model parameters are determined independently 
by fitting the measured data. However, the inferred exospheric temperature 
and absolute densities at 120 km are dependent on the existence of diffusive 
equilibrium and the assumed boundary conditions. The variations in com- 
position at 120 km are, nevertheless, a convenient way of illustrating departure 
from fixed boundary diffusive-equilibrium models; moreover, density ratios 
inferred at 120 km are less dependent on errors in exospheric and boundary 
temperatures than are the absolute densities at 120 km. 


12 



IV. MODEL PARAMETER RESULTS AND ACCURACY 


The model parameters for N 2 , O, and He densities at 450 km are given in 
Table 2 and those for O and He densities at 120 km and exospheric temperatures 
in Table 3. Also given are the estimated standard deviations for the parameters 
and a parameter indicating the overall quality of the fit: 



(8) 


where 

n‘ =i th measured density 

m 

rd = i th calculated density 

cr' 1 = estimated measurement error 

m 

N = number of data points 

d 

N p = number of parameters 

The least squares fits are carried out using about 1000 data points randomly 
selected from a data base about twice as large. 

The temperature at 120 km is fixed at 355°K and the N 2 density at 120 km 
is chosen to provide an N 2 density at 150 km very close to that recommended 
by Van Zahn [1970] and used by Jacchia [1971] , The inferred average O to N 2 
density ratio at 150 km turns out to be about 72% of that recommended by Von 
Zahn. However, in light of the errors in the various measurement methods, 
the variability of the densities, the assumptions used in extrapolating the 
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satellite measurements, and the possible error in the assumed 120 km tempera- 
ture, this small disagreement in the oxygen to nitrogen ratio is not significant. 

A total assessment of various measurements in the lower thermosphere and 
their consistency with the upper thermosphere measurements is not attempted 
in this model. 

The quality of the model fits to the data is examined by using the coefficients 
from Table 3 and plotting the measured to calculated density ratios as a function 
of time from perigee, latitude, local time, solar flux, A p , and day-count, as 
shown in Figures 3 to 5. In Figures 6 to 8 the residuals are represented by 
symbols in a latitude versus local time plot. In general, there are no large syste- 
matic differences in the fit. Examination of the plots of measured to calculated 
density ratios versus time from perigee indicates that the density data have been 
reasonably corrected for gas-surface interactions. The generally good fit does 
not necessarily mean that all significant terms have been found, particularly 
when interactions between variables are considered. For example, Figures 6 
to 8 show that while the overall distribution of points in latitude-local time space 
is good, the same is not true at all seasons. The data coverage for N 2 and O is 
much better during summer than winter and the converse is true for helium. 

This will reduce the reliability of the model in winter for N 2 and O, in the 
summer for He, and of the seasonal interaction terms. 
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The question of the validity of individual and- combinations of parameters, 
and thus the proper decomposition into the various components, is very difficult 
to answer completely. In general, one can say that since latitude profiles are 
the basic data input to the model, the model should predict latitude profiles with 
the most accuracy. While local time variations are basically derived from data 
on many different days, the lack of systematic residuals as a function of local 
time indicates that these variations should be reasonably well represented on 
average but not necessarily in their seasonal variations. The least reliable 
coefficients are expected to be those connected with the symmetrical annual and 
the symmetrical and asymmetrical semiannual components since the data base 
is not long compared to these cycles. There is also the possible confusion of 
true annual and semiannual components with these frequencies in the mean F Jo ? 
variation. While the standard deviations calculated for each coefficient give 
some guidance, they cannot completely account for a lack of uniform data 
sampling in all coordinates. 

The possibility of changes in the model coefficients from year to year was 
examined by dividing the data into two samples and fitting separately. The re- 
sults for the temperature coefficients are given in Table 4. As expected, the 
coefficients relating to the semiannual and symmetrical annual components have 
changed the most, while the majority are quite stable. 
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The probable accuracy of the inferred average temperature was tested by 
arbitrarily changing the average temperature by 50°K and 100°K and refitting 
N 2 density for the remaining coefficients. The change in Q f (Equation 8) and the 
density ratio versus time from perigee plot indicated thata50°K change would not 
be entirely unreasonable, but that a 100°K change produced fairly marked devia- 
tions from the best fit. Another consideration is that if the highest altitude den- 
sities (about three scale heights above perigee) have a systematic error of 25%, 
this would roughly translate to a 80 °K error. A test fit of oxygen density leav- 
ing T^ as a free parameter produced a temperature about 50° higher than from 
the nitrogen data. Because there is some difference between exo- 

e 

spheric temperature and the temperature near 450 km, the inferred exospheric 
temperature does depend somewhat on the s parameter. If s were .025, 
the inferred exospheric temperature would be 10° lower; and if s were .015, the 
inferred exospheric temperature would be 40° higher. The best estimate is that 
the average temperature is accurate to ±50°K. 

The accuracy with which the inferred temperature variations represent real 
temperature variations depends on the assumed constancy of the 120 km N 2 density 
and temperature. This assumption appears, however, to be a reasonable first 
approximation that can be modified later to provide better agreement with other 
measurements. A 10% variation in 120 km density would result in roughly a 
1.5% error in exospheric temperature variations and a 10% variation in 120 km 
temperature would result in roughly a 4% error in exospheric temperature 
variations. A variation of N 2 density during one day at 140 km on the order of 


16 



30% and temperature on the order of 20% was measured by Spencer et al. [1969] 
in a series of rocket flights. Semidiurnal temperature amplitudes at 120 km on 
the order of 10% have been reported by Salah and Evans [1972]. Theoretically, 
Volland and Mayr [1970] predict a 120 km diurnal temperature amplitude of 5% 
and a density amplitude of 10%. Seasonal variations in temperature at 120 km of 
about 10% have been reported by Waldteufel [1970] , with a level of relatively 
constant density at 130 km. Mayr and Volland [1972b] predict annual amplitudes 
at 120 km of 8% in temperature and 20% in N 2 density. 

The accuracy of the inferred oxygen and helium variations at 120 km depends 
on both the accuracy of the inferred temperature variations and the diffusive 
equilibrium assumption. Because of mass motions it is expected theoretically 
that there will be departures from diffusive equilibrium in the lower thermo- 
sphere, particularly for minor constituents like helium. 
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V. GLOBAL THERMOSPHERIC STRUCTURE 


The model coefficients given in Tables 2 and 3 can be used to predict at- 
mospheric structure for various conditions. Predictions for altitudes, solar 
activity conditions, etc., which are beyond the data base used in generating the 
model must obviously be used with caution. A set of plots of exospheric temp- 
erature, densities at 120 km, and densities at 450 km are presented in Figures 
9 through 14. By use of the plots for the individual components of the model a 
graphical solution for many situations can be constructed. 

The OGO-6 model shows a number of departures from previous empirical 
models based on total density measurements obtained from satellite drag [ Jacchia, 
1965; Jacchia, 1971] . Among these are the morning maximum 
for He, the high latitude of the diurnal temperature bulge position for solstice con- 
ditions, the seasonal variation of the O to N 2 density ratio with a 400°K summer 
to winter temperature difference, and the 1600 hour maximum for N 2 with an 
extended bulge in latitude. These effects will be discussed in sep- 
arate reports. The differences from drag models undoubtedly arise in part be- 
cause composition and temperature cannot be uniquely determined from total 
density measurements alone [Stein and Walker, 1965] and because of the 
smoothing inherent in the method of determining the density from variations in 
a satellite's orbit [Jacchia, 1971] . 

The total densities predicted by the OGO-6 model are compared with those 
of the J71 model [Jacchia, 1971] in Figure 15. The average density is 17% 
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higher than J71, certainly within the overall measurement error and limitations 
of model predictions. The shape of the diurnal curve for 45° latitude is in quite 
good agreement with J71 and the equatorial curve is in fair agreement, although 
the OGO-6 model has a slightly later afternoon maximum. The seasonal vari- 
ation compares well with J71, particularly north of -45° latitude. 

Measurements of temperature [Salah and Evans, 1972] made during the 
same time period as the OGO-6 measurements are compared in Figure 16. The 
average OGO-6 inferred temperature is about 7% higher than the incoherent 
scatter measurements. The annual variations are in quite good agreement and 
are over twice those predicted by J71. The diurnal variations show quite good 
agreement except for the early morning hours of winter, which is not too 
surprising considering the lack of satellite data in this area (Figure 6). 
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LOCAL TIME (hrs) 1969 I 1970 I 1971 

Figure 3 . Ratio of the measured N2 density to the model N2 density as a function of time from 
perigee, geographic latitude, local ti me, F ^ g -j solar flux, magnetic index A p and day 
of year. 
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LATITUDE (deg) 


d9e9 


LOCAL TIME (hrs) 


Ratio of the measured Nj density to the model N 2 density plotted in a latitude*loca I 
time coordinate system. The plotted letter A indicates a ratio of 1.0, B a ratio of 
l-.l, 9 a ratio of .9, and so on. In the lower half of the figure the data have been 
restricted to ±45 days from the solstices and plotted in the respective summer and 
winter hemispheres. 
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Figure 8. Ratio of the measured He density to the model He density 
plotted in the same format as Figure 6. 








Figure 9a. Contour plot of departures of inferred exospheric 
temperature from the global mean (1098 K) for Sep- 
tember equinox conditions assuming an Fjq 7 

150 and an A of 4. 
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Figure 9b. Same as 9a for June solstice conditions. 
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Figure 9c. Contour plot of departures of inferred exospheric Figure 9d. Same as 9c for June solstice conditions, 

temperature from the global mean (1098°K) using 
the diurnal, semidiurnal, and terdiurnal compo- 
nents only for September equinox conditions (as- 
suming an F-|q 7 of 150). 
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Figure 9g. Contour plot of departures of inferred exospheric 
temperature from the global mean (1098°K) using 
all annual and semiannual components (assuming 
an Fjq 7 of 150). 






Figure 10a. Contour plot of the logarithm (base e) of the ratio 
of the inferred 120 km atomic oxygen densities to 
the global mean (8.97 x 10 10 cm' 3 ) for September 
equinox conditions assuming an F 1Q 7 of 150 and 
an A of 4. 
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Figure 10b. Same as 10a for June solstice conditions. 
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Figure lOe. Contour plot of the logarithm (base e) of the ratio 
of the inferred 120 km atomic oxygen densities to 
the global mean (8.97 xlO^cm” )using the asym- 
metrical annual component only (assuming an F^g 7 
of 150). 


360 


Figure lOf. Same as lOe using the symmetrical semiannual 
component only. 





Figure lOg. Contour plot of the logarithm (base e) of the ratio 
of inferred 120 km atomic oxygen densities to the 
global mean (8.97 x 10'® cm"^) using all annual 
and semiannual components (assuming an F-|g -j 
of 150). 
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Figure lOh. Plot of the logarithm (base e) of the ratio of the 
inferred 120 km atomic oxygen densities to the 
global mean versus geographic latitude using the 
time dependent and magnetic activity components 
only. 









Figure 11c. Contour plot of the logarithm (base e) of the ratio 
of the inferred 120 km helium densities to the global 
mean (2.73 x 10 7 cm" 3 ) using the diurnal, semi- 
diurnal, and terdiurnal components only for Sep- 
tember equinox conditions (assuming an F^q 7 of 












Figure lie. Contour plot of the logarithm (base e) of the ratio 
of the inferred 120 km helium densities to the 
global mean (2.73 x 10 7 cm - ^) using the asym- 
metrical annual component only (assuming an q 7 
of 150). 



Figure Ilf. Same as lie using the symmetrical semiannual 
component only. 
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Figure llg. Contour plot of the logarithm (base e) of the ratio 
of the inferred 120 km helium densities to the 
global mean ( 2.73 x lO' 7 cm"^) using all annual 
and semiannual components (assuming an F]g 7 

of 1 50 ). 


Figure llh. Plot of the logarithm (base e) of the ratio of the 
inferred 120 km helium densities to the global 
mean versus geographic latitude using the time 
independent and magnetic activity components 
on ly. 
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Figure 12a, Contour plot of the logarithm (base e) of the ratio 
of the inferred 450 km nitrogen densities to the 
global mean (2.16 x 10^ cm~^) for September equi- 
nox conditions assuming an F^g -j of 150 and an 
A of 4, 
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Figure 12b, Same as 12a for June solstice conditions 
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Figure 12c. Contour plot of the logarithm (base e) of the ratio 
of the 450 km nitrogen densities to the global 
mean (2.16 x 1 0 6 cm - ^) using the diurnal, semi- 
diurnal, and terdiurnal components only for Sep- 
tember equinox conditions (assuming an F|q ^ of 

150). 
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Figure 12d. Same as 12c for June solstice conditions. 
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Figure 12e. Contour plot of the logarithm (base e) of the ratio 
of 450 km nitrogen densities to the global mean 
(2.16 x 1 0 6 cm -3 ) using the asymmetrical annual 
component only (assuming an F^q -j of 150). 
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Figure 12g. Contour plot of the logarithm (base e) of the ratio 
of 450 km nitrogen densities to the global mean 
(2.16 x 1 0 6 cm'A using all annual and semiannual 
components (assuming an F^q 7 of 150). 
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Figure 12h. Plot of the logarithm (base e) of the ratio of the 
450 km nitrogen densities to the global mean 
versus geographic latitude using the time inde- 
pendent and magnetic activity components only. 
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Figure 13a„ Contour plot of the logarithm (base e) of the ratio Figure 13b. Same as 13a for June solstice conaitions 

of the 450 km atomic oxygen densities to the global 
mean (6 72 v 10 7 cm -3 ) for September equinox 
conditions assuming an F,g of 150 and an A 
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Figure 13c. Contour plot of the logarithm (base e) of the ratio 
of the 450 km atomic oxygen densities to the global 
mean (6.72 x 10 7 cm' 3 ) using the diurnal, semi- 
diurnal, and terdiurnal components only for Sep- 
tember equinox conditions (assuming an Fjq 7 °f 

150). 
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Figure 13d. Same as 13c for June solstice conditions. 
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Figure 13e. Contour plot of the logarithm (base e) of the ratio 
ofthe450km atomic oxygen densities to the global 
mean (6.72 x 1 0 7 cm~^) using the asymmetrical 
annual component only (assuming an Fjq 7 of 150). 






Figure 13g. Contour plot of the logarithm (base e) of the ratio 
ofthe450km atomic oxygen densities to the global 
mean (6.72 x 10 7 cm -3 ) using all annual and semi- 
annual components (assuming an F 107 of 150). 
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Figure 13h. Plot of the logarithm (base e) of the ratio of the 
450 km atomic oxygen densities to the global 
mean versus geographic latitude using the time 
independent and magnetic activity components 
only. 
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Figure 14a.. Contour plot of the logarithm (base e) of the ratio 
of the 450 km helium densities to the global mean 
(3,02 1 0 6 cm ~ 3 ) for September equinox conditions 

assuming an F ^ q y of 150 and an A p of 4, 
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Figure 14b, Same as 14a for June solstice conditions 
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Figure 14c. Contour plot of the logarithm (base e) of the ratio Figure 14d. Same as 14c for June solstice conditions, 

of the 450 km helium densities to the global mean 
(3.02 x 10 6 cm -3 ) using the diurnal, semidiurnal, 
and terdiurnal components only for September equi- 
nox conditions (assuming an F, n 7 of 150). 
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Figure 14g. Contour plot of the logarithm (base e) of the ratio 
of the 450 km helium densities to the global mean 
(3.02 x 10 6 cm" 3 ) using all annual and semiannual 
components (assuming an F 10 j of 150). 
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Figure 14h. Plot of the logarithm (base e) of the ratio of the 
450 km helium densities to the global mean versus 
geographic latitude using the time independent 
and magnetic activity components only. 






LATITUDE (deg) 


Figure 15a. Diurnally averaged total mass density at 450 km as a function of latitude for 
solstice conditions (summer at right) for the Jacchia Model (J71) and the 
OGO-6 empirical model. 
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Figure 15b. Total mass density at 450 km as a function of local time at equinox for 
latitudes 0° and 45° for the Jacchia model (J71) and the OGO-6 empirical 
model. 
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Figure 16a. Diurnally averaged exospheric temperature as a function of 
day for incoherent scatter measurements [Salahand Evans, 
1972] and the OGO-6 empirical model. 
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Figure 16b. Diurnal variation of exospheric temperature for incoherent 
scatter measurements [Salah and Evans, 1972] and the 
OGO-6 empirical model. 



Table 1. Quiet Days with Available OGO-6 Data 



*Daycount from Jan 1 
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Table 2. Model Coefficients and Errors for 
450 km Densities 
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He 
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c 00 
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c 1 

10 

-1.551 70E+0 

3.9E-2 
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6.1 E-2 

-2.73134E-1 

1.7E-2 

O 1 
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1.8 E-2 
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10 

-5.96942E+0 
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6.15071 E-3 

3. 3 E-3 

c li 

-4.91628 E-2 
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2.79221 E-l 

3.8 E-2 
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2.5E-2 

c 1 
21 

b n 

-1.62804E-2 

1.8E-2 

3.74919E-2 

2.6 E-2 
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9.1 E-3 

-8.67463E-1 

2.5E-2 
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2.7E-2 
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1.2E-2 

b 3 1 
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-2.13267E-2 

9.7E-3 

4.75225E-2 

1.2E-2 

-9.85268E-3 

5.1 E-3 
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4.4 E-3 

-2.96619E-2 
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3.1 E-3 

d ll 
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di 
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8.2E-3 
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Table 3. Model Coefficients and Errors for Inferred 
Exospheric Temperatures and 120 km Densities 



Temperature 

He 

O 


1.09746E+3 j- 

1.1E+1 





a 20 
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n(120) 

3.0 E+ll 


2.72880E+7 

6.3E+5 

8.97076 E+ 10 

9.0E+8 

Q f 

3.05 


0.82 


1.05 


N„ 

984 


872 


905 
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Table 4. Model Coefficients for Inferred Exospheric Temperature 
Based on Different Time Periods 


n(120) 


-1.09877E-3 ± 3.5 E-l 


-1.11568E-3 ± 1.1 E-l 


n(120) 


-1 .11568 E-3 ± 1.1 E-l 


-1.09877E-3 ± 3.5E-1 



27 June 1969 

29 June 1970 


to 


to 



4 Aug 1970 

13 May 1971 

Too 

1.11568E+3 ± 

1.1E+1 

1.09877E+3 ± 

3.5E+1 

O 

o ^ 

4.55397E-3 

3.5E-3 

-8.71478 E-3 

4.1E-2 

2.59769E-2 

2.7E-3 

2.22340E-2 

1.1 E-2 

f al 
1 00 

1.12636 E-3 

4.3 E-5 

1.46310E-3 

2.3E-4 

f a 2 

1 00 

-6.14146 E-6 

1.6 E-6 

-2.34289E-5 

7.6 E-6 

f al 
00 

1.81034E-3 

4.7E-4 

2.90716E-3 

7.1 E-4 

koo 

3.67972E-3 

5.7E-4 

1.64373E-3 

6. 9 E-4 

\r a 
K 20 

4.91533E-3 

1.1 E-3 

3.79008 E-3 

1.2E-3 

coo 

2.34466 E-2 

6.0 E-3 

-1.89772E-3 

2.7E-2 

C20 

6.21293E-3 

4.5E-3 

-3.77918E-2 

3.0E-2 

*0* 

-6.97931E+1 

1.2E+1 

-2.37760E+0 

3.7E+1 

c 00 

2.65288E-2 

3.2E-3 

-1.58230 E-4 

2.1 E-2 

c 20 

1.50036 E-2 

4.3 E-3 

9.34394E-3 

3.7E-2 
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1 00 

9.69683E+1 

3.2E+0 

1.39853E+2 

1.6E+2 
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-1.63496 E-l 

5.4E-3 

-1.42320 E-l 

2.4E-2 

c 30 

-4.62570E-2 

3.8 E-3 

-5.52283E-2 

1.8E-2 

C 50 

1.68539E+2 

3.1 E-3 

1.38372E-2 

6.2E-3 


-9.81137E+0 

1.1E+0 
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3.5E+1 
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3.84413E-3 

1.0 E-2 
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2.5E-2 
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3.8E+3 
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2.8E+1 

a ll 
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a 31 
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7.8E-3 
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3.04996E-3 

6.7E-4 

2.89507E-3 

1.5 E-3 

cil 

-4.32113E-3 

7.5E-3 

8.98341 E-3 

2.0 E-2 

n 1 
C 21 

1.13675E-2 

2.1 E-3 
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8.1 E-3 

b u 

-9.75982E-2 

3.0E-3 
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3.2E-2 
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1.0 E-3 
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3.9 E-3 
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5.5E-4 
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1.3 E-3 

d ll 
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3.2E-2 
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2.0E-3 
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8. 6 E-3 
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C 

32 
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2.4E-3 

5.13250E-3 
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3.5E-3 
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4.5E-4 
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2.3E-3 
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7.9 E-4 
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1.8 E-4 
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8.6 E-4 

s 

2.05681 E-2 

9.7E-4 
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3.8E-3 

n(120) 

3.0 E+ll 


3.0 E+ll 


Q f 

2.38 


1.58 


N d 
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